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Quantized Hall conductance is a generic feature of two dimensional electronic 
systems with broken time reversal symmetry. In the quantum anomalous Hall 
state recently discovered in magnetic topological insulators, time reversal sym¬ 
metry is believed to be broken by long-range ferromagnetic order, with quan¬ 
tized resistance observed even at zero external magnetic field. Here, we use 
scanning nanoSQUID magnetic imaging to provide a direct visualization of 
the dynamics of the quantum phase transition between the two anomalous 
Hall plateaus in a Cr-doped (Bi,Sb)2Te3 thin film. Contrary to naive expecta¬ 
tions based upon macroscopic magnetometry, our measurements reveal a su¬ 
perparamagnetic state formed by weakly interacting magnetic domains with 
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a characteristic size of few tens of nanometers. The magnetic phase transi¬ 
tion occurs through random reversals of these local moments, which drive the 
electronic Hall plateau transition. Surprisingly, we find that the electronic sys¬ 
tem can in turn drive the dynamics of the magnetic system, revealing a subtle 
interplay between the two coupled quantum phase transitions. 

The integer quantum Hall effeet (QHE), first observed in elean two dimensional eleetron 
systems at high magnetie fields (1), is the paradigmatie example of a topologieal phase: differ¬ 
ent integer quantum Hall states are eharaeterized by identieal symmetries but different integer 
topologieal quantum numbers r], with the quantized Hall resistance given by R^y = rjh/e^ (2). 
However, Hall quantization may occur also in the absence of an external field as long as time- 
reversal symmetry (TRS) is broken (3). This quantum anomalous Hall (QAH) state was recently 
realized experimentally (4-9) following theoretical proposals based on combining strong spin- 
orbit coupling with long-range ferromagnetic (FM) order in magnetically doped topological 
insulators (TI) (10-18). The QH plateau transition represents the canonical example of a topo¬ 
logical phase transition, described by the divergence of the localization length and a universal 
critical scaling of transport coefficients (19). Recent theoretical calculations show that under 
certain assumptions, the QAH plateau transition can be mapped onto the same network model 
used to describe the integer QH plateau transition, leading to the same scaling laws (20). 

While the electronic topological transition in QH and QAH systems, taken in isolation, 
can be viewed as identical, the two experimental systems differ in several key aspects. The 
QH plateau transition is a purely electronic effect, in which delocalization occurs against a 
background of quenched electronic disorder. In contrast, the QAH plateau transition results 
from two coupled quantum phase transitions: the field-driven magnetic transition of the FM 
order and the electronic transition that is driven by TRS breaking of the FM transition. The 
dynamics of the FM reversal can endow the QAH electronic transition with features that do not 
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have an analog in QHE systems. Exactly at the transition, both systems can be understood as a 
network of domain walls—magnetic domain walls for the QAH, and domain walls of different 
filling factor in the QHE—that host counterpropagating chiral edge states. Crucially, however, 
in contrast to the QH plateau transition, which describes the phase diagram of the system in 
thermodynamic equilibrium, the EM domain structure in the QAH is metastable leading to 
hysteresis and relaxation dynamics that can directly affect the electronic system. Moreover, the 
energy scale of the electronic delocalization transition and its critical scaling may depend on 
the details of the microscopic EM structure and on the scaling of the magnetic phase transition. 
Einally, the electronic system, which commonly mediates EM interactions in dilute magnetic 
systems, can in turn modify the phase transition of the magnetic system. Most previous studies 
of QAH systems have used electronic transport measurements to probe the combined effect 
of magnetic and electronic evolutions, making it difficult to disentangle the individual roles of 
the two phase transitions. Exploring the reciprocal coupling of the magnetic and electronic 
quantum phase transitions thus requires selective measurement tools that can address the two 
systems independently. 

Here we combine electronic transport with a scanning superconducting quantum interfer¬ 
ence device (SQUID) of 200 nm diameter that resides on the apex of a quartz tip (SQUID- 
on-Tip (SOT)) (21, 22) to simultaneously probe the magnetic and electronic transitions at ^He 
temperatures in a 7 quintuple layer (QE) thick Cro.i(Bio, 5 Sbo. 5 )i. 9 Te 3 film grown by molecular 
beam epitaxy on a SrTiOa substrate (Eig. 1). The large dielectric constant of the substrate al¬ 
lows effective control of the chemical potential through back gating. The hysteretic transition 
between two Hall resistance (Rxy) plateaus occurs at a coercive field of He = 130 mT at 250 
mK (Eig. lA). At elevated fields, the longitudinal resistance (Rxx) shows a pronounced dip as 
a function of the back gate voltage at ~ 6 V due to an incipient QAH state (Eig. IB). At 
the lowest measurement temperature (250 mK) used here, the sample we discuss does not yet 
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show the fully developed QAH. This is similar to previous experiments on Cr-doped topolog¬ 
ical insulator films (4-9) where quantized Hall resistance appears only at dilution refrigerator 
temperatures, far below the onset temperature of hysteretic magnetic behavior (4, 23). Figures 
lE-H show images of the local distribution of the magnetic field Bz(x, y) in the sample at var¬ 
ious points along the magnetization loop. We find Bz(x, y) to be highly inhomogeneous at all 
fields, with peak contrast at He, where the average magnetization vanishes (Figs. 1F,G). Sur¬ 
prisingly, submicron structure is evident even at fields corresponding to saturation of transport 
coefficients (Figs. 1E,H). Images corresponding to opposite magnetization are highly anticor¬ 
related on microscopic scales—including at full saturation. This suggests an inhomogeneous 
distribution of magnetic moments due to segregation of the Cr dopants. While high resolution 
transmission electron microscopy measurements on our samples have yet to show any obvious 
evidence of Cr clustering (see Eig. Sll), we cannot preclude inhomogeneity at the nanoscale, 
akin to that found in other magnetically doped semiconductors such as Cr-doped ZnTe (24). 
Clear evidence of phase separation has only been seen in Cr doped Bi 2 Se 3 thin film (25) sam¬ 
ples that do not show the QAH state, while large nanoscale fluctuations in the local Cr density 
have been observed in Te based samples (26) due to random doping that introduces strong 
disorder in the material (27). 

In metallic EM thin films with out-of-plane magnetization, it is well established that magne¬ 
tization reversal develops via the nucleation and propagation of domain walls (DWs) separating 
regions of opposite magnetization. Such DW mediated magnetization reversal has also been 
imaged in ferromagnetic semiconductor films (28,29) at magnetic dopant concentrations com¬ 
parable to those in our magnetic TI films. However, scanning SOT microscopy reveals a very 
different picture of the magnetization reversal process in magnetic TIs. Eigure 2A shows a 
sequence of Bz(x, y) images acquired for increasing values of yoH near He. The five images 
appear almost identical; however, numerical subtraction of successive image data (ABz(x, y). 
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see Fig. 2B) reveals the underlying dynamic process. Instead of the anticipated DW motion, 
magnetization reversal occurs through a series of random events in which isolated nanoscale is¬ 
lands undergo a reversal of their out-of-plane magnetic moment (see Movie SI). As we discuss 
below, this constitutes a direct microscopic observation of superparamagnetism in magnetically 
doped TI films. Our observations caution against drawing conclusions about the ferromagnetic 
state solely from macroscopic magnetization probes (SQUID magnetometry, magneto-optical 
Kerr effect) which show square hysteresis loops with robust zero field remanence (30). 

To quantify the superparamagnetic dynamics across the Hall plateau transition, we fit each 
of the local features in the ABz(x, y) maps with a point-like out-of-plane magnetic moment m 
(Fig. S3). Figure 2C summarizes ~1700 such fits, accumulated over four different ranges of 
magnetic field. Throughout the measured range, the spatial distribution of reversal events is ran¬ 
dom (Fig. 2D inset and Fig. S4), suggesting weak interactions between neighboring islands and 
supporting the conclusion that superparamagnetism in Cr-doped BiSbTe films is characterized 
by aggregations of dopant atoms. The deviation of the Hall resistance from its saturation value 
begins around yoH = 0 and is accompanied by the reversal of nanoscale moments with average 
ffi = 3 X 10^/iB (Figs. 2C,D). Given an average saturation magnetization of ~ 3y^/Cx atom as 
obtained from global magnetization measurements (Fig. S14), the estimated average diameter 
of these flipping islands is d = 51 nm for our 7 nm thick film, considerably below our spatial 
resolution of ~300 nm. As the field is increased towards He, however, a pronounced change in 
the moment distribution is observed with a shift to higher m values and appearance of a long 
tail of large moments with m > IO^/ib (Figs. 2C,D). The microscopic reversal moments m can 
be summed to obtain the net change in magnetization M over a continuous field range (Fig. 
2E). Comparison with simultaneously acquired R^y shows a qualitative match, implying that 
the behavior of the transport coefficients through the plateau transition is mainly determined by 
the underlying change in magnetization. Superparamagnetic behavior and a similar relationship 
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between measured magnetization and transport coefficients were found in a second sample as 
well as in a Mn-doped Bi 2 Te 3 film (see Movie S2 and Figs. S2, S12-13 and S15). 

Hysteretic magnetic transitions are a signature of metastability, and typically display tem¬ 
poral relaxation via thermal activation or quantum tunneling. We probe magnetic relaxation in 
real time by polarizing the system at hqH = — 1 T and then ramping the field to a positive set 
point fioHset- We then acquire repeated images of Bz{x,y) while simultaneously monitoring 
electronic transport. No relaxation is observed on laboratory time scales for yoHset < 0 in 
either magnetization or transport coefficients. Spontaneous relaxation begins to be evident at 
small positive fields (yoHset = 63 mT), manifesting as magnetic reversal events ABz{x, y) and 
a slow upward drift of R^y. The frequency and number of these reversals increases significantly 
near the coercive field at yoHset = 126 mT (Figs. 3A,C). 

The temporal relaxation measurements further corroborate the superparamagnetic behavior: 
at temperatures well below the blocking temperature, the magnetization of a superparamag- 
net is hysteretic, showing minimal relaxation at low fields. On approaching He the magnetic 
anisotropy barrier U is reduced, leading to relaxation when U ~ k^T. Since U is proportional 
to the volume of the superparamagnetic particles, smaller islands undergo thermal activation 
at a lower field. Simultaneous transport measurements (Fig. 3D) indicate that the electronic 
transition closely tracks the magnetic relaxation, with transport coefficient relaxation evident 
at yoHset = 63 mT and pronounced at 126 mT, in accord with the total temporal change in 
magnetization extracted from the SOT data (Fig. 3C). 

The plateau transition observed in electronic transport appears to be mainly dictated by the 
underlying magnetic reversal. Surprisingly, however, we find that the dynamics of the magnetic 
system can in turn be influenced by the electronic system. To explore this effect, we perform a 
sequence of magnetic imaging at 126 mT interspersing consecutive scans with small excursions 
of the back gate voltage AV^ (Fig. 3B). Remarkably, even small AVg ~ 1 V excursions enhance 
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the relaxation of the superparamagnetic islands significantly, as is evident from the statistics of 
the observed moment reversals (Fig. 3C). This enhancement in turn is evident in the transport 
coefficients, as shown in Fig. 3D. 

A full comparison of the effects of applied magnetic field, gate voltage, and time on transport 
coefficients is presented in Fig. 4C, which shows a unified parametric plot of the transport 
coefficient vector {Rxx, Rxy)- On this plot, the large magnetic hysteresis evident in Figs. 4A,B 
is absent, demonstrating that the relation between Rxx and Rxy (at a given Vg) is a universal 
function determined by the magnetization. Within a single constant-arc-shaped plot, zero 
net magnetization corresponds to the maximum of the arc at zero Hall angle, while the varying 
Hall angle along the arc reflects the varying sample magnetization. Variable Vg traces over 
AVg = ±30 V at hqH = ±1 T, marked in black, show contours of variable carrier density at 
constant saturated magnetization. Note that at full saturation, variable Vg data retraces the same 
path upon repetition. 

On this plot, the temporal relaxation of transport coefficients (Fig. S8) is seen to track 
constant Vg arcs, implying that temporal magnetic relaxation is the dominant mechanism and 
consistent with the scanning magnetometry results (Figs. 4C,D, gray dots along = 6 V arc). 
Again consistent with magnetometry data, Vg excursions enhance the magnetization relaxation 
dramatically in the metastable regime (Figs. 4 C,D). The blue trace shows the evolution of 
the transport coefficients during three Vg excursions with AV^ increasing from ±1 V to ±3 
V at a constant field of 126 mT. Successive Vg sweeps do not retrace each other, consistent 
with an irreversible change in the net magnetization. Repeating the same experiment for larger 
excursion of AV^ = ±30 V, the magnetization is observed to rapidly relax, nearly reaching 
complete saturation (cyan in Figs. 4C,D). As Fig. 4C makes clear, the gate-induced magnetic 
relaxation is strongly dependent on maximum extent of the voltage excursion, although we find 
it to be independent of sweep direction and rate. 
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Strikingly, despite the dramatic effect of gate voltage variations on magnetic relaxation, field 
sweeps at different values of Vg show only small deviations in the coercive field (Figs. 4A,B) 
indicating that carrier density has little effect on the average magnetization. These seemingly 
contradictory observations can be qualitatively understood by invoking the strongly disordered 
nature of the superparamagnetic state. In our strongly disordered system, the global gap in the 
density of states is bridged by a proliferation of subgap states induced by the random polar¬ 
ization of the superparamagnetic islands. These localized states can in turn mediate the FM, 
in addition to the proposed global van Vleck mechanism that arises from band states (15). Vg 
excursions may modify the magnetic anisotropy energy of individual islands through the strong 
dependence of local density of states on position and energy, randomly changing the magnetic 
potential landscape without significantly changing the average magnetization and coercive field. 
Unlike in QHE systems, however, this disorder is not quenched: the weakly interacting super¬ 
paramagnetic islands are metastable. At any positive field, some local configurations are sep¬ 
arated from the ground state by small energy barriers U, and any local perturbation can cause 
irreversible flipping of some of the islands leading to rapid stimulated relaxation upon repeated 
excursions of Vg. This mechanism is of course enhanced near H^. Our statistical analysis of the 
flipping moments at various Vg in Fig. S5 qualitatively supports this mechanism. 

The observed superparamagnetic state and its rich dynamics emphasize the intricate cou¬ 
pling between the electronic and magnetic sub-systems in magnetic TIs. Our results suggest 
that the quantum phase transition between QAH states can be strongly affected by the nature of 
the underlying magnetic transition leading to deviations from the expected universal scaling of 
the plateau transition as indicated by recent transport studies {31). Even though the topological 
state is robust to the presence of magnetic disorder at very low temperatures, nanoscale super¬ 
paramagnetism may be responsible for the fragility of the QAH state at elevated temperatures. 
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Figure 1: Electrical transport and scanning magnetic imaging of 7 QL thick 

Cro,i(Bio.5Sbo.5)i,9Te3 sample at T =250 mK. (A-B) Transport measurements showing mag- 
netie field dependenee of R^x (red) and Rxy (blaek) at = 6 V (A) and the Vg dependenee at 1 
T (B). The dip in Rxx marked by the arrow shows the ineipient QAH state. (C) An optieal image 
of the sample and SOT showing the eleetrieal contacts and the SOT reflection from the sample 
surface. (D) Electron micrograph of the SOT used for the magnetic imaging. (E-H) Scanning 
SOT images 5x5 /im^ of the out-of-plane magnetic field Bz{x, y) at ~ 300 nm above the sam¬ 
ple surface at four anti-symmetric locations along the magnetization loop marked in (A). Note 
strong anti-correlation between (E) and (H), and (F) and (G). Pixel size 50 nm, pixel integration 
time 10 ms. 
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Figure 2: Magnetization reversal dynamics. (A) Sequence of SOT magnetic images Bz{x, y) 
taken at consecutive magnetic fields in 0.5 mT steps at T = 250 mK. (B) Differential images 
ABz{x, y) obtained by subtracting pairs of consecutive Bz{x, y) images in (A) showing the iso¬ 
lated magnetic reversal events (red) of the superparamagnetic moments (see Movie SI). (C) Sta¬ 
tistical analysis of 1690 reversal events attained over ranges of magnetic fields centered around 
four hqH values: total number of moment reversals Nm, average magnetic moment m, average 
superparamagnetic island diameter d, and rate of the magnetization change dM/d^yoH) over 
the given range. (D) Chart of relative contribution of different moment sizes m to the total 
magnetization change M within two field ranges centered at = —15 mT (yellow) and 
HqH = 154 mT (blue). Inset: location of the moment reversals within the field range around 
HqH = 154 mT. (E) Cumulative magnetization change M due to moment reversals m over four 
field ranges (left axis, colored symbols) and the simultaneously acquired R^y (right axis, black). 
The total magnetization in each range is offset by an arbitrary constant. 
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Figure 3: Temporal and back gate induced relaxation of the superparamagnetic state. (A) 

Differential image ABz{x,y) obtained by subtraetion of two eonseeutive images aequired at 
eonstant /loHset = 126 mT and = 6 V after a field ramp from - IT. Image acquisition time 
is 200 sec with 50 sec wait time between images. (B) Same as (A) with gate excursion progres¬ 
sively increasing from AV^ = 0.1 to 1.1 V in-between consecutive images. (C) Histogram of 
the temporal relaxation process showing the moment reversals m attained from four consecu¬ 
tive ABz{x,y) images at y-oHset = 63 mT (dark blue) and at yoHset = 126 mT (light blue), 
and of V^-induced relaxation at y^Hget = 126 mT acquired following the temporal relaxation 
of 20 minutes (green). (D) R^y as a function of field (black) and during relaxation at fixed field 
taken simultaneously with the magnetic imaging. Temporal relaxation over 20 minutes is more 
pronounced at 126 mT (light blue) than at 63 mT (dark blue). Vg excursions (green) induce 
large relaxation at 126 mT. Inset: full Rxy hysteresis loop showing the region of interest. 
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Figure 4: Transport measurements and universal plot of magnetic relaxation. (A) R^y and 

(B) vs. applied field at T = 250 mK and different Vg showing magnetie hysteresis with 
similar Hf.. (C) The same data plotted as universal are-like eurves of R^x vs. Rxy at various 
Vg. Extrema of the arcs correspond to saturation magnetization at -1 T (-1-1 T) on the lower left 
(right) end of each arc. Gray dots indicate 60 min of temporal relaxation at /loHset = 126 mT 
and Eg = 6 V (see also Fig. S9). Gate sweeps at ^qH = ±1 T (black lines) trace the ends of 
the arcs, and are reversible. Gate sweep at 126 mT (blue and cyan) are metastable, inducing 
magnetic relaxation and propagation along the arcs from Rxy < 0 towards positive saturation. 
(D) The Rxy relaxation data (gray, blue, cyan) and the Rxy field sweep at Eg = 6 V (green). 
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1 Movies of the magnetic moment reversal dynamics 


The fabrication of the SOT devices and the scanning SOT magnetic imaging were performed 
as described in Refs. 21,22. The scanning magnetic images of Bz{x, y) and the differences be¬ 
tween consecutive images ABz{x, y) were acquired at constant increments of applied magnetic 
field and compiled into movies (Movies SI and S2). Figures SI and S2 show representative 
frames from the movies of two Cr-doped (Bi,Sb) 2 Te 3 samples. 
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Figure SI: Movie snapshot of magnetization reversal process in 7QL Cro i(Bio. 5 Sbo. 5 )i. 9 Te 3 
film at T=250 mK. (A) Scanning SOT image 5x5 of the out-of-plane magnetic field 
Bz{x,y) at ~ 350 nm above the sample surface. The running title indicates the value of the 
applied field which is increased in steps of 0.5 mT between consecutive images (see Movie SI). 
Image acquisition time was 200 seconds with pixel size of 50 nm. The median background 
value of the signal is subtracted in each image. (B) Differential image y) obtained by 

subtracting the preceding y) image from the current one in (A). Isolated magnetic reversal 
events of the superparamagnetic moments are visible as red-green peaks. The characteristic size 
of these peaks is determined by the size of our SOT (200 nm) and by the height of the SOT above 
the sample, while the physical size of the flipping FM islands is only few tens of nm (see Fig. 
2) and is below our resolution. The moments flip in the positive direction of the ascending 
applied field giving rise to the positive value of the local ABz{x, y) field peaks. Note the much 
smaller span of the color code in (B) as compared to (A). (C) R^y of the sample measured 
simultaneously with the magnetic imaging. The red circle shows the current value while the 
black circles are the data at previous field values. (D) Same as (C) for R^^:- Figures 2A and 2B 
in the main text show few images from the movie while the statistical analysis of the flipping 
moments in the movie are presented in Figs. 2C-E for field range /iqH = 154 mT. 
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Figure S2: Movie snapshot of magnetization reversal process in lOQL Cr-doped 
(Bi,Sb) 2 Te 3 film at T=250 mK. (A) Scanning SOT image 5x5 of the out-of-plane mag¬ 
netic field Bz{x,y) similar to Fig. SI but on a different sample and on descending magnetic 
field in steps of 1 mT (see Movie S2). Image acquisition time was 180 sec with pixel size of 50 
nm. (B) Differential image ABz{x, y) obtained by subtracting a previous Bz{x, y) image from 
the one in (A). The superparamagnetic moments flip in the direction of the descending magnetic 
field giving rise to the negative values of the local peaks in AB^ix, y). (C) R^y of the sample 
measured simultaneously with the magnetic imaging. (D) Same as (C) for R^x- 
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2 Magnetic moment fitting procedure 


The differential images ABz{x, y) were fitted by a distribution of point-like magnetie moments 
with out-of-plane polarization as deseribed in Fig. S3. 




Numerical fit 
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Figure S3: Fitting procedure of magnetic moments. (A) Differential image ABz{x, y) in 7QL 
Cro.i(Bio, 5 Sbo. 5 )i. 9 Te 3 film obtained by subtraeting Bz{x,y) image at = 154.0 mT from 
image at yoH = 154.5 mT (same as the bottom image in Fig. 2B). Numerieal proeedure is then 
used to find loeal peaks marked by x that are above a threshold of 4 yT set by the measurement 
noise level. (B) Best fit proeedure is then applied in whieh eaeh peak (five in this ease) is fitted 
to the magnetie field generated by a point dipole 2m at height h below the SOT eonvoluted with 
the aetive area of the SOT (the faetor 2 arises from the fact that the superparamagnetic island 
magnetization flips from —m to +m). (C) The quality of the fit is demonstrated by subtracting 
the fit (B) from the image (A). (D) The resulting magnitudes m of the fitted dipoles at a fixed 
optimal h = 345 nm. The diameter d of the superparamagnetic islands is then calculated from 
their moment m assuming an average saturation magnetization of Sys/Cr atom (obtained from 
global magnetization measurements presented in Fig. S14), Cr doping of 5%, unit cell size of 
0.43 nm, and thickness of 7 nm. The resulting size of the islands is smaller than SOT diameter 
or the height h justifying the point dipole aproximation. A full simulation taking into account 
the finite size of the islands gave similar results. 
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3 Spatial distribution of the magnetization reversals 


MqH = -15mT 



MqH = 80mT 



IJ^H = 133mT 



IJqH = 154mT 



Figure S4: Spatial distribution of the magnetization reversal process in 7QL 

Cro.i(Bio. 5 Sbo. 5 )i. 9 Te 3 sample at T=250 mK. (A-D) Location of all the moments that flipped 
upon increasing the applied field from ^qH = —30 mT to 1 mT (A), ^qH = 60 to 100 mT 
(B), hqH = 126 to 140 mT (C), and /iqH = 140 to 167 mT (D) corresponding to the four field 
ranges in Figs. 2C-E. The random distribution of the moment locations further demonstrates 
that the magnetization transition occurs through reversal of weakly interacting superparamag- 
netic islands rather than through FM domain wall motion. 
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4 Vg dependence of the dynamics of magnetic relaxation 

Figures 4A,B show that Vg has only a small effect on the coercive field He while AVg excursions 
near He relax the magnetization dramatically as shown in Figs. 4C,D. We have also investigated 
the temporal relaxation process at different constant values of Vg as shown in Fig. S5. 
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Figure S5: Statistical analysis of temporal moment relaxation for different Vg values in 

7QL Cro.i(Bio sSbo sli.gTes film at T=250 mK. The histogram shows the relative contribution 
of moments of different magnitudes m to the total magnetization relaxation M for three values 
of Vg. For each histogram the magnetic state was initialized at a fixed Vg by sweeping the 
applied field from -1 T to 125 mT and consecutive y) images were acquired for 25 min 
at fixed fioHset = 125 mT similar to Fig. 3A. While for Vg = —23 and 6 V the distributions 
are similar, at = 37 V a clear shift of the distribution to larger moments is visible showing 
that at the same applied field the dynamics of the relaxation process has some dependence on 
Vg. These results indicate that the magnetic anisotropy energy of individual islands is affected 
by Vg apparently through variations in the disorder-induced random local potential or density 
of states. Since the local variations are random their average effect on the global coercive field 
could be small as observed in Figs. 4A,B. 
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5 Transport measurements 


The samples were mechanically patterned into a Hall-bar geometry and the transport measure¬ 
ments were performed using three synchronized SR830 lock-in amplifiers for simultaneous 
measurement of the current, longitudinal voltage, and transverse voltage (Fig. S6). The applied 
current was 10 to 20 nA at a frequency of 17.71 Hz. In order to minimize the mixing between 
longitudinal and transverse resistances switching between the current and voltage leads in VdP 
configuration was performed using Keithly7001 switch box for most of the presented data. The 
gate voltage was applied using a DC source (Keithley2400). The sweeping rate of Vg was lim¬ 
ited to 0.5 V/sec to ensure absence of heating, and after reaching the target Vg there was a 3 
seconds waiting time before the measurement. 

The simultaneous magnetotransport and SOT imaging measurements were performed in 
Oxford Heliox ^He cryostat equipped with 12 T magnet. In order to prevent heating the mag¬ 
netic field ramping was kept at a low constant rate of 0.4 T/min in all experiments. This rate 
was used for magnetotransport measurements as well as for preparation of the initial magneti¬ 
zation state for temporal and gate-induced magnetic relaxation measurements. Figures S7 to S9 
show the dependence of transport coefficients on temperature, back gate voltage, and the tem¬ 
poral relaxation. Since large AVg excursions cause pronounced relaxation of the magnetization 
we have also compared the hysteretic behavior of the transport coefficients in presence of gate 
excursions as shown in Fig. SIO. 
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Figure S6: Schematics of transport measurements. Three panels showing the seheme of the 
seratehed Hall bar appearing in Fig. 1C. The width of the eentral ehannel is 500 /im, the inner 
distanee between two voltage probes on the same side is 1000 /im and the width of eaeh voltage 
ehannel is 200 fim. The orange part marks the two eurrent leads and the ehannel between them 
i.e. where bulk eurrent flows. The red lines indicate measurement configuration of Vxx or 
Vxy. (A) Hall bar measurement configuration for simultaneous measurements of I, Vxx and Vxy 
using three lock-in amplifiers resulting in Rxx = ^ and Rxy = (B) VdP configuration for 

measurement of Rxx with switching between current and voltage leads using a switch box. (C) 
VdP configuration with switching for measuring Rxy. 
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Figure S7: Temperature dependence of the transport coefficients in 7 QL 

Cro.i(Bio.5Sbo.5)i.9Te3 film. Hysteretic loops of R^x and Rxy vs. applied field at 1/^ = 8 V at 
two temperatures T = 240 mK and 1.5 K. The eoereive field at 1.5 K is redueed to about a 
third of its value at the base temperature. 
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Figure S 8 : Gate voltage dependence of R^x in 7 QL Cro i(Bio sSbo.s)! 9 X 03 sample. The 

longitudinal resistance Rxx at iiqH = 1 T vs. 1^ at T = 250 mK and 1.6 K. The incipient QAH 
state is clearly visible at T = 250 mK as a pronounced dip around Vg = 8 Y marked by the 
arrow. A smaller dip is resolved also at T = 1.6 K. 
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Figure S9: Temporal relaxation of transport coefficients near the coercive field in 7 QL 

Cro.i(Bio.5Sbo.5)i.9Te3 sample. Time dependenee of and R^y at HoHset = 126 mT at 
= 6 V after ramping the field from -IT. The relation between R^^ and R^y follows the 
universal are-like sealing as shown in Fig. 4C by the gray dots. 
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Figure SIO: Transport coefficients in 7QL Cro.i(Bio sSbo 5 )i. 9 Te 3 sample with continuous 
gate excursions. Hysteretic loops of and R^y vs. applied field at 1/^ = 6 V without (blue) 
and with (red) gate exeursions. In the former ease the transport eoeffieients were measured 
eontinuously while sweeping the applied field. In the latter case at each value of the applied 
field, gate excursion of Al/^ = ±30 V was applied for ten times prior to the measurement of 
R^x and R^y Increasing the number of AV^ repetitions did not induce additional change in the 
transport coefficients. The effective coercive field is seen to be reduced substantially by the gate 
excursions in contrast to Fig. 4A,B which shows no appreciable change in He for field sweeps 
at various constant gate voltages over the same range of Vg values. Note that the resulting 
reduction in is much larger than what can be achieved by temporal relaxation (see Figs. S9 
and 4D). Despite the large reduction in the hysteresis by the gate excursions a fully reversible 
equilibrium state could not be reached. 
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6 Growth and structural characterization 


Samples were grown by moleeular beam epitaxy (MBE) on eommereial SrTiOa substrates with 
dimensions of 5 x 5 x 0.5 mm^. Substrates were annealed in oxygen at 925° C for 2 hour 
30 min typically and were screened using atomic force microscopy (AFM) to make sure the 
surface was atomically ordered with an RMS roughness of less than 0.3 nm. Substrates were 
indium mounted for growth in the MBE chamber and outgassed at a substrate temperature of 
550° C for about 1 hour to remove any residual contamination before growth. The TI films were 
grown using elemental materials of at least 5N purity evaporated from Knudsen type thermal 
cells at a growth rate of about 0.4 QE per minute at a temperature of 250° C, as measured by a 
pyrometer. 

Transmission electron microscopy (TEM) studies were carried out on a sample grown under 
similar conditions. The TEM sample was prepared by focused ion beam (EIB) and measured 
in an EEI Titan double aberation corrected TEM at 200 keV. High angle annular dark field 
scanning TEM (HAADE) images showed a well ordered crystal with the expected quintuple 
layer structure. Two slightly bright lines in the middle of each QE correspond to the location 
of the heavier Bi atoms. A disordered amorphous region about 1 nm thick was observed at the 
interface with the SrTiOa. We have observed a similar layer in undoped films as well. Energy 
dispersive spectroscopy (EDS) was used to map the spatial distribution of elements. Within 
each of six regions studied the distributions of elements appeared uniform with the exception 
that Bi appeared deficient in the interfacial region. 
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Figure Sll: STEM imaging and EDS elemental mapping of a Cr doped (Bi,Sb) 2 Te 3 film 
on SrTiOs. (A) HAADF image showing the erystal strueture and EDS elemental maps of the 
various atomie speeies. Less Bi is observed in the disordered layer at the interfaee with the 
SrTiOs- (B) Linecut made by averaging over the area shown. The Te-Bi(Sb)-Te-Bi(Sb)-Te 
QL structure is clearly seen. Te, Bi, Sb and Cr are shown in green, cyan, magenta, and blue 
respectively. (C) Large scale HAADL image of the film. 
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7 Transport and superparamagnetic dynamics in lOQL Cr- 
doped [Bi, Sb) 2 Tes sample 


Similar studies were performed on an additional lOQL Cr-doped sample. Figures S2, S12, and 


S13 show that the transport eoeffieients in this sample are far from QAH quantization, however 


the superparamagnetie behavior is qualitatively the same. 



Figure S12: Magnetization reversal dynamics and transport coefficients in lOQL thick 
Cr-doped (Bi,Sb) 2 Te 3 thin film at T=250 mK. (A) Sequence of 5 x 5 /rm^ scanning SOT 
magnetic images Bz{x,y) taken at consecutive descending magnetic fields in 1 mT steps. (B) 
Differential images ABz{x, y) obtained by subtracting pairs of consecutive Bz{x, y) images in 
(A) showing the isolated negative magnetic reversal events of the superparamagnetic moments 
(blue). A larger sequence of images is available in Movie S2 as described in Fig. S2. (C) 
Transport coefficients vs. applied field at three values of Vg. 
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Figure S13: Scaling of the cumulative magnetization change and the transverse resistance 
in lOQL Cr-doped (Bi,Sh) 2 Te 3 film at T=250 mK. Evolution of magnetization M (left axis, 
blue cireles) on deseending magnetie field (after sweep down from 1 T) obtained by eumulative 
addition of the superparamagnetie moment reversals M = 2m, derived from a sequenee of 

images like in Fig. S12 and the simultaneously measured Rxy (right axis, red). The data points 
were aequired only at ranges of fields where the SOT had high sensitivity resulting in the gaps 
between the ranges. The magnetization in eaeh range is offset by an arbitrary eonstant. Through 
most of the transition region a elose relation between R^y and M is observed. 


18 



8 Global magnetization studies 


Magnetization measurements were taken of a 40 nm thiek sample using a eommereial SQUID 
magnetometer with the applied field perpendieular to the film (Fig. S14). For the temperature 
sweeps the sample was cooled from 305K to 3.5K in either a 9000 Oe field for the field cooled 
(FC) curve or 1.3 Oe field for the zero field cooled (ZFC) curve (1.3 Oe is the remenance in 
the magnet after demagnetization, so this corresponds to zero field). The applied field is then 
set to 50 Oe and the magnetization is measured as the sample is warmed. The same procedure 
was followed using an annealed STO substrate with no film. The measurements from the STO 
substrate were subtracted from the measurements of the film to remove the contribution from 
the substrate. After this a small constant offset was subtracted from the FC and ZFC curves so 
that the magnetization is zero at 50K. A clear peak is seen in the ZFC data below T^ at around 7 
K, suggestive of a blocking temperature associated with superparamagnetism. This peak is also 
seen in thin (~10 QL thick) samples as well as in samples grown on InP(l 11)A substrates. Such 
peaks are typically associated with superparamagnetic or spin glass materials but have also been 
occasionally observed in ferromagnets. Further studies are needed to identify its origin. 

For the field sweep, the sample is cooled in field to 5K and the applied field is swept from 
9000 Oe to -9000 Oe and back. Again, the substrate contribution is subtracted. All measured 
magnetizations are normalized by the sample area. It should be noted that the features in the 
field sweep at -I-/-2500 Oe are artifacts due to the small magnetization signal in the raw data 
when the strong diamagnetic contribution of the substrate causes the total magnetization to cross 
through zero at these values. Using a composition of Cro.i(Bio, 5 Sbo, 5 )i. 9 Te 3 and a high field 
saturation magnetization of 7.25e-7 emu/mm^ gives a moment of approximately 3.15 /ig/Cr, 
close the expected value of 3. 

Given these results, a 5x5 iirn? area of a 7QL sample of Cro.i(Bio. 5 Sbo. 5 )i. 9 Te 3 should show 
a magnetization change of ~3.3 x when integrating from zero magnetization to saturation 
value. This number is comparable to M obtained by cumulative addition of m in Fig. 2E. For 
lOQL sample the total change in M from negative to positive saturation should be ~9.4 x10Vb 
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comparable to Fig. S13. 
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Figure S14: SQUID magnetometry measurements of a ~40QL thick film. (A) M vs. T field 
cool and zero field cool magnetization curves showing a peak in the ZFC below Tc- (B) M vs. 
H hystersis curve at 5 K after subtracting the STO substrate eontribution. 
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9 Transport and superparamagnetism in Mn-doped BiTe 
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Figure S15: Transport and scanning magnetic imaging of 70 nm thick Mn-doped BiTe 
sample at T=250 mK. (A) R^y vs. applied field showing hysteretie loop with eoereive field 
of He — 0.2 T in 70 nm thiek Mn-doped BiTe sample grown by MBE on InP substrate. (B- 
F,I) Seanning SOT images 2.5 x 5 /rm^ of the magnetie field Bz{x, y) at six points along the 
magnetization loop marked in (A). Images (B) and (C) were aequired near zero applied field (45 
mT and 23.5 mT respeetively) on the opposite branehes of the loop. The magnetie strueture is 
highly antieorrelated similar to the behavior in Fig. 1. Conseeutive images at fields of 180 (D), 
195 (E), and 200 mT (E) show superparamagnetie behavior with no domain wall motion. (G) 
Differential image ABz{x,y) obtained by subtraeting (D) from (E) showing isolated moment 
reversals. (H) Same as (G) for images (E) and (E). (I) Bz{x, y) at 0.6 T - nearly at the closing of 
the loop - demonstrating that though with lower contrast, the magnetic nano-structure prevails 
up to high fields. 22 









10 Atomic Force Microscopy of Cr-doped samples 


10QL Cr-doped sample 7QL Cr-doped sample 



Figure S16: AFM topography images of the two studied Cr-doped (Bi,Sh) 2 Te 3 samples. 

Both samples show surfaee disorder with eharaeteristie lateral length seale of about 100 nm. 
X-ray diffraetion shows that this surfaee roughness is ereated by the formation of a mosaie of 
twinned erystalline domains during epitaxial growth. The 7QL sample deseribed in the main 
text (right eolumn) has a more uniform morphology and lower surfaee roughness than the lOQL 
sample. 
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